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AERONAUTIC SYMBOLS 

1. FUNDAMENTAL AND DERIVED UNITS 



Length 
Time__ 
Force. _ 

Power _ 
Speed _ 



Symbol 



I 
t 

F 



P 
V 



Metric 



Unit 



meter 

second 

weight of 1 kilogram. 



horsepower (metric) . 
fkilometers per hour. 
\meters per second, _ 



Abbrevia- 
tion 



m 
s 

kg 



k.p.h. 
m.p.s. 



English 



Unit 



foot (or mile) 

second (or hour)__ 
weight of 1 pound 

horsepower 

miles per hour 

feet per second 



Abbrevia- 
tion 



ft. (or mi.) 
sec. (or hr.) 
lb. 



hp. 

m.p.h. 
f.p.s. 



2. GENERAL SYMBOLS 



9> 

m, 

i 



Weigh t=mg 

Standard acceleration of gravity =9.80665 
m/s 2 or 32.1740 ft./sec. 2 
W 

Mass== — 

Moment of inertia —mk 2 . (Indicate axis of 

radius of gyration k by proper subscript.) 
Coefficient of viscosity 



p, Kinematic viscosity 
p, Density (mass per unit volume) 
Standard density of dry air, 0.12497 



kg-m" 4 -s 2 at 



^ 15° C. and 760 mm; or 0.002378 lb.-ft." 4 sec. 2 
Specific weight of "standard" air, 1.2255 kg/m 3 
0.07651 lb./cu. ft, 



or 



3. AERODYNAMIC SYMBOLS 



b 

c, 
b 2 
S' 
V, 

2, 

L, 

D, 

Do, 

D<, 

D p , 

C, 

B, 



Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio 

True air speed 

Dynamic pressure=^p 1 P r2 

Lift, absolute coefficient C L =^ 



Drag, absolute coefficient C D = 



D 



Do 



Profile drag, absolute coefficient C^o^^g 
Induced drag, absolute coefficient C Di =^ 
Parasite drag, absolute coefficient ^> P = ^~g 

C 

Cross-wind force, absolute coefficient C c =^-^ 
Resultant force 



Q, 

VI 



#0> 



7, 



Angle of setting of wings (relative to thrust 
line) 

Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 

Reynolds Number, where I is a linear dimension 
(e.g., for a model airfoil 3 in. chord, 100 
m.p.h. normal pressure at 15° C, the cor- 
responding number is 234,000; or for a model 
of 10 cm chord, 40 m.p.s., the corresponding 
number is 274,000) 

Center-of-pressure coefficient (ratio of distance 
of c.p. from leading edge to chord length) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 
Flight-path angle 
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DETERMINING PROFILE DRAG 

By Harry J. Goett 



SUMMARY 

An experimental investigation has been conducted in the 
full-scale tunnel to determine the accuracy of the Jones and 
the Betz equations for computing profile drag from total- 
and static-pressure surveys in the wake of wings. Sur- 
veys were made behind 6- by 36-foot airfoils of the N. A. 
C. A. 0009, 0012, and 0018 sections at zero lift and behind 
the N. A. C. A. 0012 airfoil at positive lifts. The surveys 
were made at various spanwise positions and at distances 
behind the airfoil ranging from 0.05c to 3.00c. 

The reduction of the test data by either the Jones or the 
Betz equation gave profile-drag coefficients agreeing within 
2 percent with those obtained by force tests at zero lift. 
The variation of the profile drag determined at stations 
from 0.05c to 3.00c behind the trailing edge was small and 
the error resulting from the induced field of a lifting airfoil 
did not exceed 2.5 percent at a C L of 1.0 and a spanwise 
station of 0.78 b/2. 

INTRODUCTION 

The use of the momentum met hod for the detennina- 
tion of profile drag has recently increased, owing- mainly 
to the equations developed by Betz (reference 1) and 
by Jones (reference 2) by which the method has been 
made applicable in the region of increased static pres- 
sure close behind a body. The derivation of these 
equations, which arc based on the original principle 
stated by Froude in 1874, requires certain assumptions. 
The errors introduced by these assumptions have been 
the subject of theoretical analysts (references 2 and 3), 
which have set an upper limit for the ('n'(n'< involved 
but fail to define their actual value. 

The investigation reported herein was conducted to 
determine experimentally the magnitude of these errors 
by determining I he el Feci of a number of variables upon 
the measured drag. (See reference 4.) The necessary 
wake surveys were made in the N. A. C. A. full-scale 
wind tunnel behind symmetrical airfoils of three thick- 
ness ratios. The effect of distance behind the airfoils 
was first investigated by a comparison of drag deter- 
minations made at locations ranging from 0.05c to 
3.00c behind (he trailing edge. A check was then 

obtained on the accuracy of the method by a comparisoD 



with force-test dra^ measurements at zero lift. Fin- 
ally, the effect of the induced-flow system of a lifting 
wing was investigated. 

SYMBOLS 



The symbols used in the report are defined as follows: 
Hq, free-stream total pressure. 

total pressures in field of airfoil. 

(See fig. 4.) 
free-stream static pressure, 
static pressures in field of airfoil, 
free-stream dynamic pressure, l/2pt V- 
free-stream velocity, 
local velocity in field of airfoil, 
hypothetical velocity in wake (Betz 

equation), 
vertical coordinate of point, 
airfoil chord, 
elemental areas, 
density, 
airfoil span. 

velocity along the Y axis, 
velocity along the Z axis, 
airfoil profile drag, 
airfoil profile-dra^ coedicient . 
section profile-drag coefficient, 
airfoil lift coefficient, 
sect ion lift coefficient. 



fli, //- Hz 
Vo 

Pl> ft, P3 

<Zo 

Uo. 

Uo' 

V 
c 

dS, dSt, dS 2 , dS 3 

p 
b 

v, 
w 



@D 0 
ci 



EXPERIMENTAL INVESTIGATION 

APPARATUS 

The experimental work was conducted in the N. A. 
C. A. full-scale wind tunnel (reference 5). This tunnel 
has a turbulence factor of U as determined by sphere 
tests (reference 6). A typical static-pressure gradient 
along the axis of the tunnel (jet empty) is shown in 
figure 1. This gradient was allowed for in determining 
the free-stream reference pressure for the momentum 
measurements. The buoyancy effect of the gradient 
is small. 

Three 6- by 36-foot rectangular airfoils having 
N. A. C. A. 0009, 0012, and 0018 sections were used in 
these tests. The airfoils, which were covered with 
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K6-inch aluminum slieets, had all the screw heads filled 
and the surface painted, sanded, and polished to a 
glossy waxlike finish to insure aerodynamic smooth- 
ness. The airfoil tips were rounded, each tip forming 
one-half of a solid of revolution with the radius at each 
chordwise station equal to one-half the local airfoil 
thickness. Figure 2 shows one of the airfoils mounted 
in the t unnel jet. 
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Figure l- 



-2-10 12 3 

Distance from trailing edge, chords 

-Typical static-pressure gradient along test section of full-scale wind 
tunnel (jet empty). 



The rack used for the total- and the static-pressure 
surveys consisted of a comb of 39 total-pressure tubes 
and one of 13 static-pressure tubes. These combs 
were spaced G inches laterally and the entire assembly 
Was mounted on the survey carriage. The detailed 
spacing and the dimensions of the tubes on both combs 
are shown in figure 3. Each tube was connected to a 
multiple-tube, photographic-recording manometer car- 
ried in the survey carriage. 




f f 




FIGURE 2.— The 6- by 30-foot N. A. C. A. 0012 airfoil mounted in the tull-scule wind 

i unnel. 



A survey was first made with the total- and the static- 
pressure combs at each station of measurement with 
the jet empty. This survey established the total- 
pressure and the static-pressure gradients in the 
tunnel at the points of measurement. Pitch-angle 
surveys were next made behind the airfoils to establish 
the average down wash angle across the held of measure- 
ment. Total-pressure and static-pressure readings were 
then taken in the wake region with the rack perpendicu- 
lar (o the average downwash direction at each station. 



This procedure kept the effect of flow angularity on the 
measurements at a minimum, since the local angle 
across the rack varied no more than ±3° from the 
average. The effect of periodic pressure fluctuations 
in the tunnel jet was eliminated by the instantaneous 
readings taken on the photographic manometer. 

SCOPE OF TESTS 

Pressure and drag measurements were made at loca- 
tions and under conditions as follows: 

I. Total- and static-pressure surveys were made at 
zero lift behind the three airfoils at 27 span wise loca- 



To multiple - 
tube 

manometer 



\ I 



To multiple- 
tube 

manometer 



4" 



6" 



24" 2\m_, 



i 

2-i 



6' 



.065' O.D. x .036" ID. 
Comb of total - 
pressure tubes 



Four '/&"D.. 
static - 



-/" 



pressure , 
orifices"'' .125" ID. 
Comb of static 
pressure tubes 



Fig u he 3.— Combs of total- and static-pressure tubes. 

tions, 0.15c behind the trailing edge. At the 0.06 6/2 
station, surveys were obtained at longitudinal stations 
varying from 0.05c to 3.00c behind the trailing edge. 
Force tests were made to furnish comparative drag 
data. 

2. Total- and static-pressure surveys were obtained 
behind the N. A. C. A. 0012 airfoil at lift coefficients of 
0, 0.28, 0.47, 0.65, 0.83, and 1.13 at six spanwise loca- 
tions, 0.15c and 0.30c behind the trailing edge. Force 
tests were made to furnish comparative drag data. 

All tests were run at an air speed of 90 miles per hour, 
giving a test Reynolds Number of 5,000,000. 

THEORY 

The profile drag of a body can be determined from 
the loss of momentum per unit time that it imposes upon 
the free stream. If a region exists behind the body 
where the static pressure has returned to that of the 
free stream (fig. 4 (c)), the profile drag of a nonlifting 
body will be given by the expression 

Do= Pf w f U z (Uo- UJdS ( 1 ) 

where IT indicates that the integration is confined to 
the wake region. For practical reasons, it is desirable 
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in most cases to make the survey in the region close 
behind the airfoil where p is in excess of p$ (fig. 4 (b)). 
In this region the drag will be equal to 

D Q =ff(ih + P mdSt - ff (p 2 + P W)dS 2 (2) 

where bolli integrations are carried to infinity. 

Since it is impossible to survey to infinity as required 
by equation (2), this equation must be transformed into 
one involving only quantities in the wake region. This 
transformation has been made by Betz (reference 1) 
and by Jones (reference 2). 



Pi HrH> 



H 0 JJo ,p 0 



(a) 




Figure 4.— Diagram of airfoil and wake. 



Betz builds up a hypothetical How by means of a 
system of sources of such strength that the total 
pressure in the wake of the body is restored to the 
value it would have in potential (low. (See fig. 4, 
H 2 =H 0 , U 2 =U 2 .) This system has a resultant 
thrust equal to the thrust of the sources. It differs 
from the real system only in the region of the wake so 
that the difference in thrust between the two systems 
is equal to the difference in momentum per unit time 
passing through the wake region of each. It then 
follows that 

Drag of real system = (Difference in thrust between 
hypothetical and real systems) — (Thrust of hypo- 
thetical system) 

Thus, the integration over the region external to the 

wake is eliminated and the expression for profile drag 

reduces to 

D 0 = J'" f(H Q -H 2 )dS 

+| ({U 2 '-U 2 ){U 2 '+U 2 -2U 0 )dS (3) 

In terms of total and static pressures to be measured, 
the section profile-drag coefficient becomes 



i mrHo-Hf] r 



H 0 — p 0 



Jl 



(4) 



The last step involves the assumption that, for a lifting 
airfoil, 

II 2 =\pU 2 *+p 2 (5) 

instead of 

Ih=\p{U 2 2 +v 2 2 +w 2 >) +p 2 (6) 

Jones (reference 2) assumes a hypothetical flow in 
the wake in which there is no energy interchange 
between tubes of How in the wake behind the point of 
measurement and consequently that Bernoulli's equa- 
tion may be applied to these, tubes of (low. On the 
basis of such an assumption, the total and the static 
pressures measured close behind a body in a region of 
increased static pressure give sufficient data to obtain 
the corresponding velocity loss (and therefore mo- 
mentum defect) at a point where the static pressure 
has reached the free-stream value. Tims on a non- 
lifting airfoil, 

r> 0 =pf"fu 2 (u«-Ui)ds (7) 

and (on the foregoing assumption) 

H2-Po=IpU 3 2 (8) 

Then 

2/ "f^H^2^W I Vo- jH^jftdS (9) 
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Reduced to coefficient form, equation (9) becomes 

CJ ^H 0 —po\ ^Ho — p () / 

which Jones also applies to a lifting airfoil. 

The effect of the assumptions made in the derivation 
of the Betz and the Jones equations has received con- 
siderable study. The errors involved in the method of 
Betz are difficult to estimate and the validity of the 
derivation is difficult to establish. 

Taylor (reference 3) has shown thai the neglect in 
the Jones method of tin 4 internal tangential stresses 
("mixing") which occur in the wake downstream of the 
measured section is theoretically unsound. From the 
examination of a number of typical profiles, Taylor has 
shown that the error does not exceed 1.5 percent but he 
also shows that much larger errors are possible. 

The induced field of a finite lifting wing may cause 
errors in the methods. First, the assumption made in 
equation (5), that the v and the w components may be 
neglected, will be a source of error. Second, there is 
the possibility that the vortices in the wake region may 
damp out causing a loss of total pressure, which appears 
erroneously as profile drag. An analysis of this possi- 
bility, based on certain typical wake profiles, has been 
made by Jones in reference 2 and the maximum value of 
the error due to this possible pressure loss has been esti- 
mated. 
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In addition to the foregoing errors, inaccuracies will 
possibly arise from incorrect readings of static and 
total prrssuros caused by turbulence and stream angu- 
larities behind the airfoil. 

RESULTS AND DISCUSSION 

COMPUTATION OF RESULTS 

In the computation of results, the values of II 2 and 
p 2 across the wake profile were determined from faired 
curves of total and static pressures, to which a correc- 
tion was applied to allow for the vertical gradients 
existing in the tunnel. The values of 77 0 and j>o were 
determined from readings taken well outside the w ake 
with a correction applied to obtain the values of these 
quantities at the position of the airfoil. These values 
were then substituted in equations (4) and (10) and the 
results were plotted against the vertical position in the 
wake. The resulting curve was integrated, the sum- 
mation being the section profile-drag coefficient at the 
station of measurement. An additional correction was 
applied for displacement of the effective center of the 
total-pressure tubes in a velocity gradient. 

EFFECT OF DISTANCE BEHIND THE AIRFOIL 

The variation in measured drag with distance behind 
the airfoil is shown in figure 5 for the three airfoils at a 
c t of 0.05 and for the N. A. C. A. 0012 at values of 
Ci of 0.78 and 1.32. Each point is the average of results 
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Figure 5.— Variation of measured drag with distance behind the airfoil. 

from two or more surveys. The curves show that 
the greatest difference in drag, as measured at the fore- 
most and the rearmost positions, is approximately 3 
percent. This difference is within the experimental 
scatter of the measurements, estimated to vary from 
±1 percent at the 0.15c station to ±3 percent at the 
3.00c station (where the wake profiles are shallow and 



wide). It is therefore concluded that the measured 
drag, as evaluated hy either the Jones or the Betz 
method, is unaffected by distance behind the airfoil 
within the accuracy of the measurements. 

Figure 5 also indicates that there is no significant 
difference between the drag as determined by the Betz 
and the Jones equations. The maximum spread 
between the two methods is less than 1 percent. All 
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Figure 6.— Results of drag determinations obtained by the momentum method by 
Muttray, Kramer and Doetsch, and Jones. 

further drag determinations were therefore made by the 
Jones equation because of the greater simplicity of the 
necessary computations. 

The results of other similar investigations (references 
2, 7, and 8) arc shown in figure 6. From these results, 
Muttray concluded (reference 7) that sufficient data 
had not yet been obtained to warrant the conclusion 
that the measured drag was independent of the distance 
behind the airfoil. The present investigation is con- 
sidered, however, to have furnished sufficient data to 
support this conclusion. 

EFFECT OF TURBULENCE ON STATIC-PRESSURE MEASUREMENTS 

In relation to the possible effect of turbulence on 
the measurement of static pressure, in the wake, a 
comparison of measured static pressures with computed 
static pressures behind the N. A. C. A. 0012 and 0018 
airfoils at zero lift is shown in figure 7. The pressures 
behind the airfoils at zero lift were computed for the 
case of ideal flow about the airfoils and for flow with 
a boundary Layer and a wake by means of a source-sink 
distribution to represent the airfoil and the wake. 

COMPARISON OF MOMENTUM- AND FORCE-TFST RESULTS 

The accuracy of the momentum method is indicated 
from a direct comparison with force- test results. The 
drag coefficients obtained from momentum surveys at 
27 span wise locations at 0.15c behind the three airfoils 
at zero lift are plotted in figure 8. These curves, when 
integrated across the span, give an over-all C Dq for 
each airfoil. The drag coefficients obtained in this 
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TABLE I 

COMPARISON OF PROFILE-DRAG COEFFICIENTS AT 



ZERO LIFT OBTAINED 
FORCE TESTS 



FROM MOMENTUM AND 





C/> 0 
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FlGUBE 7.— Experimental and theoretical static-pressure variation behind the N. A. C. A. 0012 
and 0018 airfoils at zero iift. 



.0120 



O 
o 

V 



o/oo 



0080 



0060 



^.0040 
0. 



0020 



CO 



















1 1 1 

A/ A r. A 






























O 




- 0009 
■ 00/2 - 
















r 












A 


00/8 


















































X- 




I 

\ 












r — 










N 




* 1 
















e 








~~ -x- 






c— -yi 






— X— 


— x-~ 


— J-* 

V \ 








„-^> — 




n 


-i 


>— U 





















































































































































































































cient against section lift coefficient obtained under the 
various conditions of induced flow that exist between 
t ho center line and the tip of a finite airfoil at positive 
lifts. Such curves are given in figure 9 for six spanwise 
locations from 0.06 to 0.90 b/2. The section lift coef- 
ficients were computed on the basis of a lift distribution 
given by Glauert (reference 9), which was found to check 
well with pressure-distribution tests. These curves, 
having been slightly shifted in order to make them 
agree at zero lift, are superimposed in figure 10. It 
will be noted that, out to 0.78 b/2 and up to a c x of 1.0, 
the maximum dispersion is ±2.5 percent from 
the 0.06 6/2 curve. This variation compares 
with a difference of approximately 7 percent 
that is indicated by Jones' analysis, which is 
based on the assumption of complete damping 
out of the vortices (reference 2). Inasmuch 
as part of the dispersion in figure 10 is due 
to experimental scatter, the ±2.5 percent is 
considered a conservative estimate of the 
effect of the induced field upon the measure- 
ments. At the 0.90 b/2 station, the distorted 
curve indicates that, above a section lift co- 
efficient of zero, the air-stream angularities 
become such as to make the measurements 
unreliable. 
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Figuke 8.— Variation of section profile-drag coefficient across the span of the rounded-tip airfoils 

at zero lift. 



manner are compared with those measured by force 
tests in table I. The maximum difference is less than 
2 percent, indicating the order of accuracy of the 
momentum method. A similar comparison cannot be 
made at positive lifts because it is impossible to obtain 
results from the momentum method in the region of 
the airfoil tips owing to the intensity of the vortices. 

EFFECT OF INDUCED FLOW 

The total effect of the errors caused by the induced- 
flow field of a lifting airfoil may be determined from a 
comparison of the curves of section profile-drag coeffi- 



A direct comparison between drag results 
obtained from momentum and force tests of 
a lifting airfoil is impossible because of the 
failure of the momentum method near the 
tip as well as the inclusion of induced drag 
in the force-test measurements. If the in- 
duced drag is deducted from the force-test 
drag, however, the two methods should give 
results differing only by the drag contributed 
by the tips. Such a comparison has been 
made in figures 1 1 and 12. A plot of profile 
drag (determined by the momentum method ) 
against span-wise position is given in figure 
1 1 ; the curves were extrapolated in the tip 
region and no allowance was made for an in- 
crease in drag at the tips. Integration of these 
curves across the span gives the average pro- 
file-drag coefficient C Dq , which has been plotted against 
C L in figure 12. The result is compared with the profile- 
drag coefficient determined from force tests in the usual 
manner (i. e., by deducting the computed induced drag). 
A curve of section characteristics obtained by the 
momentum method (0.06 b/2 curve from fig. 10) is also 
given to show the comparison with the average profile 
drag across the span. 

An appreciable spread will be noted between the 
momentum- and the force-test results; the difference 
varies from 1 percent at zero lift to 22 percent at a 
C L of 1.0. A number of causes other than tip effects 
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may contribute to this difference, including possible 
errors in the momentum method and uncertainties in 
computed induced drag and in the several corrections 
applied to the force-test results. 

CONCLUSIONS 

The results presented herein lead to the following 
conclusions with regard to the determination of profile 
drag by the momentum method and the application of 
the Betz and the Jones equations under the conditions 
of the present investigation: 

1. The drag determined by the momentum method 
did not vary appreciably with distance behind the air- 
foil between stations ranging from 0.05c to 3.00c 
behind the trailing edge. 

2. At zero lift, the drag determined by the mo- 
mentum method agreed with that measured by force 
tests within 2 percent. 

3. Inboard of 78 percent of the semispan, the effects 
of the induced-flow system of a lifting wing did not 
cause errors exceeding 2.5 percent at a c t of 1.0. 

4. The Betz and the Jones equations gave results 
that agree within 0.5 percent at stations ranging from 
0.05c to 3.00c back of the trailing edge. 

5. For measurements made no farther than 3.00c 
behind the trailing edge, the experimental scatter 
varied from 1 percent at zero lift to 3 percent at a 
d of 1.0. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., December 20, 1938. 
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Absolute coefficients of moment 

r - A r -M. 

(rolling) (pitching) 



(yawing) 



Angle of set of control surface (relative to neutral 
position), 6. (Indicate surface by proper subscript.) 



4. PROPELLER SYMBOLS 



D, 

P> 

P/D, 

V, 

V„ 

T, 

Q, 



Diameter 
Geometric pitch 
Pitch ratio 
Inflow velocity 
Slipstream velocity 

T 

Thrust, absolute coefficient (7 T = ^ 2 p 4 

' Q 



Power, absolute coefficient Cp= ^ 3 ^5 

Speed-power coefficient^-^ pj^ 
Efficiency 

Revolutions per second, r.p.s. 
Effective helix angle =tan'^2~j 



Torque, absolute coefficient C Q = 



5. NUMERICAL RELATIONS 



1 hp.=76.04 kg-m/s=550 ft-lb./sec. 
1 metric horsepower= 1.0132 hp. 
1 m.p.h. =0.4470 m.p.s. 
1 m.p.s.=2.2369 m.p.h. 



1 lb. =0.4536 kg. 

1 kg=2.2046 lb. 

1 mi. = l,609.35 m=5,280 ft. 

1 m=3.2808 ft. 



